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A nonequilibrium molecular dynamics computer simulation on microsegregated solutions of sym- 
metrical diblock copolymers is reported. As the polymer concentration increases, the system under- 
goes phase transitions in the following order: body centered cubic (BCC) micelles, hexagonal (HEX) 
cylinders, gyroid (GYR) bicontinuous networks, and lamellae (L), which are the same morphologies 
that have been reported for block copolymer melts. Structural classification is based on the patterns 
of the anisotropic static structure factor and characteristic 3-dimensional images. The systems in the 
BCC micellar (per 3 = 0.3) and HEX cylindrical (pa 3 = 0.4) phases were then subjected to a steady 
planar shear flow. In weak shear flow, the segregated domains in both systems tend to rearrange 
into sliding parallel close-packed layers with their normal in the direction of the shear gradient. At 
higher shear rates both systems adopt a perpendicular lamellar structure with the normal along 
the neutral direction. A further increase in the shear rate results in a decrease in lamellar spacing 
without any further structural transitions. Two critical shear rate values that correspond to the 
demarcation of different structural behaviors were found. 

I. INTRODUCTION 

Block copolymers are macromolecules composed of se- 
quences, or blocks, of chemically distinct repeat units. 
Due to the low entropy of macromolecules, even a small 
distinction between the units causes, below the order- 
disorder transition (ODT) temperature, a phase sep- 
aration. Because of the chemical connectivity of the 
blocks, block copolymers cannot separate into macro- 
scopic phases. Instead, spatially periodic nanostructures 
are formed. The structure and its symmetry, or morphol- 
ogy, is determined by the volume fraction of blocks, or 
ultimately by the copolymer composition. 

We shall consider the simplest type of block copoly- 
mers: linear AB diblock copolymers. In this case, 
the copolymer composition is characterized by the ra- 
tio / of the number of monomeric units in the block of 
species A, Na, to the overall degree of polymerization, 
N = N A + N B : f = N A /N. We would like to character- 
ize the possible equilibrium morphologies of block copoly- 
mer systems by referring to block copolymer melts, which 
have been fairly well studied both experimentally and 
theoretically. 1 ' 2 Sorted in ascending /, these morpholo- 
gies are cubic (BCC or FCC) packing of spherical micelles 
(S), hexagonally packed cylinders (C), bicontinuous cu- 
bic ("gyroid") phase (G) , and lamellae (L), with the last 
corresponding to the symmetrical case of / = 1/2. 

When diblock copolymers are dissolved in a selective 
solvent, i.e. a solvent favorable for one component (e.g., 
A) and poor for the other (B), the morphologies also 
depend on the concentration of the block copolymers. 
The full phase diagram has been successfully obtained 
experimentally 3 and confirmed by self-consistent mean- 
field (SCMF) theory. 4 Interestingly, the morphological 
changes that occur by varying the concentration of di- 
block copolymers in a selective solution are essentially 
the same as those that occur in melt by varying the 



composition. Above the critical micelle concentration 
(CMC), symmetrical diblock copolymers in a selective 
solvent form spherical micelles which, with an increase 
in concentration, arrange into cubic (BCC or FCC) mi- 
cellar crystals, and then successively transform to hexag- 
onally ordered cylinders, bicontinuous networks, and, fi- 
nally, lamellae (Disordered micelles — » S — » C — > G — » 
L) . The similarity of the two scenarios can be understood 
by the following simple consideration. The addition of a 
selective solvent (or increasing the selectivity of the sol- 
vent) to the melt of block copolymers corresponds to an 
increase in the volume fraction of the selected block, and 
thus the phase behavior of the solution as a function of 
the concentration can be mapped onto the melt with in- 
creasingly asymmetric copolymer compositions. 

Effects of shear on block copolymer systems include 
shear-induced ordering, domain reodering, and disorder- 
ing. In the isotropic-lamellar transition, by suppressing 
nonlinear fluctuations, shear flow raises the transition 
temperature; thus, in a certain temperature region the 
lamellar phase can be induced by applying shear. The 
lamellar ordering occurs with wave vector normal to both 
shear and velocity gradient. 5 ' 6 

Reordering of block copolymer nanostructures to new 
low-frictional phases occurs when the flow becomes faster 
than the relaxation rate of the system, for a review of re- 
cent experimental reports see Ref. [7,8]. The lamellar 
phase under a weak planar shear flow (low shear rate 
or low shearing frequency) adopts an orientation paral- 
lel to the shear plane, with the lamellar normal along the 
flow gradient direction, whereas at higher shear rates, the 
lamellae flip to the perpendicular orientation, with their 
normal along the flow vorticity (neutral) direction. 9 A 
further increase in the shear rate results in a reduction in 
lamellar spacing. When the hexagonal cylindrical phase 
is subjected to shear flow, the cylinders are perfectly ori- 
ented along the flow direction, with the hexagonal plane 
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{10} either parallel, at low shear rates, or perpendicular, 
at higher shear rates and near ODT, to the flow. It has 
also been shown that spatial fluctuations of the cylindri- 
cal domains, if characterized by a lifetime greater than 
the inverse shear rate, lead to break-up, or disordering, 
which has in fact been observed at higher shear rates. 
Cubic micellar phases under shear flow change to more 
complicated symmetries, such as twinned BCC and FCC 
lattices and hexagonal close-packed (HCP) layers, with 
their densest planes (e.g. {110} for BCC and {111} for 
FCC) slipping in parallel to the shear plane. At higher 
shear rates, a loss of long-range order due to the shear- 
induced micelle melting has been reported. 7 ' 8 

Although the reorganization and reorientation of do- 
mains in microphase-segregated block copolymer systems 
have been reported, there have been few studies on the 
breakup of domains. Since domains remain intact and 
due to the aforementioned mapping, previous results re- 
garding block copolymers under shear can be reviewed 
without separating the results related to melts or sol- 
vents. We sought to study block copolymers in a selective 
solvent under a high enough shear so that we could ob- 
serve domain breakup and the corresponding new phase 
transitions. 

A mesoscopic simulation technique based on the dy- 
namic mean-field density functional theory has been re- 
cently employed to investigate lamellar orientation in 
block copolymers subjected to shear and was able to cap- 
ture parallel and perpendicular lamellar states at low and 
higher shear rates. 10 

In this work, we use the non-equilibrium molecular 
dynamics (NEMD) 11,12 computer simulation method to 
model block copolymers under shear flow conditions and 
investigate structural changes in the systems. The data 
regarding rheological and other properties obtained for 
the same systems will be published elsewhere. 



II. MODEL AND METHOD 

A. Model potentials 

Using the bead-spring model, we studied systems of n 
symmetrical diblock copolymer chains with N = Na + 
Nb beads per chain and a composition / = Na/N = 0.5. 
The selectivity of the solvent is modeled using different 
parameters for the Lennard- Jones (LJ) potential. The 
bead pairs A — A and A — B interact via purely repul- 
sive, truncated at the minimum and shifted LJ potential 
(Weeks-Chandler- Andersen potential) : 
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Since e can be used to set the energy scale and a can be 
used for the length scale, we will henceforth use so-called 
Lennard-Jones reduced units where e = a = 1 and also 
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FIG. 1: Geometry and kinematics of the planar shear flow 



set the mass m of the particles to unity, so that time 
is measured in units of (cr 2 m/e) 1 / 2 . B — B interaction 
includes long-ranged attraction and is modeled by the 
following form of LJ potential, modified so as to satisfy 
the continuity and smoothness requirements at the trun- 
cation point: 
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The radius of interaction is chosen to be sufficiently small, 
r c = 2, to avoid possible "freezing" of multi-particle ag- 
gregates and to increase the speed of calculation. The 
potential has a minimum at r « 2 1 / 6 ; a = 4.913708 gives 
4> as the well depth. We take the temperature and well 
depth as T = 1 and <j>/T = 1.2, which implies that beads 
A and B are in favorable and unfavorable solvent condi- 
tions, respectively. 

For any two adjacent beads along the chain, in addi- 
tion to the excluded volume LJ potential, Eq. (1), the at- 
tractive anharmonic FENE (finitely extensible non-linear 
elastic) spring potential is also used with spring constant 
k = 30 and maximum extension Rq — 1.5: 
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B. Equations of motion 

We model the systems under conditions of steady 
planar shear flow (viscometric flow) 13 . The geometry 
and kinetics of the flow are illustrated in Fig. 1. The 
phase space is sampled by integrating the thermostat- 
ted SLLOD equations of motion 11 for each particle i G 
[l,nN]: 



qi = Pi/m + i x jq yi , 
pi = Fi - i x jp yi - Xpi 

dx (t) = 7, 



(4) 
(5) 
(6) 
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where the position of the particle i is qi(q X i,q y i,q Z i) T , 
the momentum is 'Pi{PxiiPyiiPzi) T i Fj is the force act- 
ing on the i th particle, and d x is the lattice strain as- 
sociated with the LeesEdwards periodic boundary con- 
ditions. The parameter A is the Gaussian thermostat 
multiplier: 
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The Gaussian thermostat fixes the instantaneous kinetic 
energy along a trajectory, K(t) = l/2m^™=iPi • Pi = 
K(0) = K . If the popular Gear predictor-corrector 
method 14 ' 15 is used to numerically integrate such equa- 
tions of motion, the kinetic energy tends to drift away 
from its initial value. Therefore, a newer symplectic- 
like integrator obtained in Ref. [16] using an operator- 
splitting technique is employed, with the integration time 
step as large as 0.01, which is one order of magnitude 
greater than in Gear method. 



C. Technical details 

Simulations were performed in a cubic cell of size 
L = 40 with periodic boundary conditions for systems 
of symmetrical diblock copolymers with a fixed chain 
length of N = 10+10. Systems with different concen- 
trations in the range pa 3 = 0.2 ~ 0.8 were simulated to 
obtain a zero-shear phase diagram. The effect of shear 
flow on systems with two concentrations was considered: 
per 3 = nN/L 3 = 0.3 with the number of chains n — 960, 
and pa 3 = 0.4 with n = 1280. Without flow, the first 
block copolymer system corresponds to a BCC lattice 
of spherical micelles (S), and the second to hexagonally 
packed cylinders (C). For each shear rate, simulations 
were carried out from a random configuration until the 
time dependencies of energy and other variables of in- 
terest appeared to be constant. The typical times to re- 
quired to reach a steady state varied from as long as 
to = 6 x 10 4 at low shear rates to to = 4 x 10 3 at the 
highest shear rate 7 = 2. 

Solvent molecules surrounding block copolymers are 
not simulated explicitly in our model. Instead, the main 
phenomcnological effects of the medium such as shear 
impulse transfer and thermostating arc included in the 
equations of motion, and selectivity is ensured by the 
effective potentials. Such a mean-field approach to the 
medium corresponds to the situation when the time and 
space scales of processes in the solvent are much smaller 
than those of the solute. Thus, the values calculated 
represent the corresponding polymer contributions. 



TABLE I: Phases of symmetrical diblock copolymers in a se- 
lective solvent. DS - disordered spherical micelles; S - BCC 
lattice of spherical micelles; C - hexagonally arranged cylin- 
ders; G - bicontinuous, "gyroid" , phase; L - lamellar phase. 
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where the summation is performed over Ns B-beads that 
are inside the circle inscribed in the simulation cell. Thus, 
Sb corresponds to a scattering experiment that probes 
only the order of the B-rich domain structure. The an- 
gular brackets denote an ensemble average that also cm- 
braces the non-constant Ns- For the anisotropic systems 
studied, the scattering patterns depend naturally on the 
direction of the beam. By choosing the scattering wave 
vector q appropriately, the patterns on the three orthog- 
onal planes, with normal along the shear, shear gradient, 
and neutral directions, were obtained. 



III. RESULTS 

A. Zero-shear phases 

The morphological transitions caused by varying the 
concentration of symmetrical diblock copolymers in se- 
lective solvent in the absence of shear flow are schemati- 
cally given in Table I. The diagram qualitatively agrees 
with the previously published experimental 3 , theoretical 4 
and lattice Monte Carlo computer simulation 18 ' 19 results. 
With an increase in concentration, the system adopts the 
following sequence of morphologies: DS— ^S^C^G 
— ► L. Therefore, the model and method used in this work 
can be concluded to be suitable for studying structural 
changes in block copolymers. However, some expertise 
and special techniques, e.g. simulating annealing when 
the temperature of the system being equilibrated is grad- 
ually reduced, are required to attain equilibrium for the 
systems without shear flow. These difficulties can be at- 
tributed to the large number of particles (up to as many 
as 51200), the associative character of interactions, and 
the effects of the periodic boundary conditions. In con- 
trast, as will be described later, a steady state is attained 
much faster and easier under shear flow. 



D. Static structure factor 



Effect of shear flow 



The quantity that helps to examine the structure 
is the static structure factor, which gives the position 



The structural changes caused by shear flow are illus- 
trated in Figs. 2 and 3. Each row in both figures contains 
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FIG. 3: Typical snapshots and scattering patterns for the system pa 3 — 0.4. 
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a typical real-space snapshot of the configuration of B- 
rich domains, with A-rich domains in the interspacings 
omitted for clarity. The next three subfigurcs in a row 
are gray-scale maps of the intensity of the structure factor 
(scattering patterns) projected on the (zy) (shear direc- 
tion), (xz) (shear gradient) and (xy) (neutral direction) 
planes. Different rows correspond to the steady state 
structures observed for different shear rates. They show 
the evolution, from the top to bottom, of the systems 
with increasing shear rates. 



1. BCC micellar crystal sheared 

The results for a polymer density of 0.3 are presented in 
Fig. 2. At zero shear rate the system exhibits a BCC lat- 
tice of spherical micelles. With the application of shear, 
the simple cubic symmetry is broken and six-fold patterns 
appear on the zy pattern (b). Such a six- fold pattern has 
been observed experimentally 20 and was explained as lo- 
cal distortion of the (110) BCC plane into 2-d HCP lay- 
ers. The layers also reorient, with their normal perpen- 
dicular to the shear gradient (c). Thus, the flow proceeds 
via parallel sliding of the closely packed layers. With an 
increase in the shear rate above 7* w 0.05, the sliding 
parallel layer structure becomes unstable and the phase 
transition commences (d), leading to the formation of 
perpendicular lamellae. The new structure can be distin- 
guished unambiguously at 7** w 0.4, and further changes 
include stabilization and better separation within lammc- 
lae (e) and a shear-induced decrease in lamellar spacing 
(f)- 

2. HEX cylinders sheared 

The results for a density of 0.4 are presented in Fiq. 3. 
While it is difficult to achieve a cylindrical phase in HEX 
ordering in the absence of shear flow, even a small shear 
leads to a much faster formation of HEX cylinders aligned 
along the shear direction, (a). The structure remains sta- 
ble after shear flow is turned off and differs from that ob- 
tained without any shear only in terms of the orientation 
of cylinders. The energy of the state is also lower than 
that of any apparently metastable state generated with- 
out pre-shearing. This suggests that such pre-shearing 
techniques can be used to facilitate attaining equilibrium 
structures in complex fluids like microsegregated block 
copolymer systems. 

Greater shear breaks the percolating cylinders into 
rod-like micelles. An increase in the shear rate up to 
the first critical value 7* causes the system to rearrange 
into sliding parallel layers composed of the rod-like mi- 
celles (b). An increase in the shear rate above 7* causes 
the onset of the transition (c) to perpendicular lamellae, 
which is completed by the second critical shear rate 7** 
(d). A stronger shear leads to enhancement of the nor- 
mal density fluctuations of lamellae and their splitting, 
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FIG. 4: Typical snapshots of chain conformations in a perpen- 
dicular lamellar structure. The chains lie in the (xz) plane, 
and therefore it is sufficient to only show this projection. 



or to shear-induced decrease in lamellar spacing, which 
is apparent in the snapshot of the simulation cell as the 
splitting of four into five lamellae. 



C. Chain conformations 

Computer simulation of polymeric systems makes it 
possible to directly observe chain conformations. In the 
systems studied, due to the selective solvent, the con- 
formations of blocks of A and B species differ. In the 
absence of shear flow, ^4-blocks are in swollen coil-like 
conformations while £?-blocks are in a rather compact 
globule- like state. Shear flow is known to stimulates poly- 
mer chains to elongate and orient in the shear direction 21 . 
We observed such a behavior for both ^4-blocks and B- 
blocks, with flow overcoming the associative character 
of the interaction of the latter. In view of the forma- 
tion of perpendicular lamellae under high flow, where 
block copolymer conformations have to be compatible 
with lamellae spacing in z direction, a question arises of 
whether the chain orientation is still along the direction 
of shear. Typical chain conformations in the perpen- 
dicular lamellar phase are shown in Fig. 4 in the (xz) 
plane. This projection shows essentially all of the in- 
formation needed because the chains on average lie in 
the shear plane. The figure shows that both subchains 
flip between two directions on average oriented along the 
shear direction. This observation is supported by studies 
of the subchain gyration tensor which will be published 
elsewhere. 



IV. CONCLUSIONS 

The non-equilibrium molecular dynamics method was 
applied to study block copolymer systems in a selec- 
tive solvent under a shear flow field. The systems un- 
dergo structural changes that were studied with regard 
to the patterns of the asymmetric static structure factor 
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TABLE II: Structural behaviour of diblock copolymers in se- 
lective solvents under shear flow depending on the shear rate. 
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and typical 3-dimensional snapshots. The structure fac- 
tor, unlike the instantaneous snapshots, provides a more 
reliable basis for drawing conclusions, and showing the 
patterns in all three mutually perpendicular projections 
makes it possible to classify the new morphologies. The 
results enable us to clearly distinguish two trends, as 
shown in Table II. In weak flow, the systems show par- 
allel sliding of close-packed layers, with the normal along 
the shear gradient. At higher shear rates, the more fa- 
vorable morphology is perpendicular lamellae, with the 
normal along the neutral direction. This behavior is ob- 
served for both systems with two different polymer con- 
centrations. 

The following explanation can be adopted from studies 
of block copolymer melts, where the parallel and perpen- 
dicular orientation of the lamellar phase under weak and 
strong shear flow, respectively, has been investigated 9,22 . 
In the parallel alignment, less viscous strata permit the 
easy sliding of mechanically contrasting layers, which in- 
creases the basin of attraction of this steady-state ori- 



entation. However, near the ODT, where the mechan- 
ical contrast is not so large, stronger shear leads to 
perpendicular orientation. It has been argued that a 
steady flow suppresses concentration fluctuations and a 
weak shear more strongly suppresses fluctuations along 
the shear gradient, thus stabilizing parallel orientation. 
While a stronger shear more strongly suppresses fluctu- 
ations along the vorticity axis. 

In our study, we observed that flow via parallel and 
perpendicularly oriented layers or lamellae is also seen 
for block copolymer solutions. 

Two critical values of the shear rate are found. One, 
7* « 0.05, corresponds to the end of the region of the 
stability of parallel sliding layers and the onset of the 
transition leading to perpendicular lamellae, which arc 
clearly formed at 7** « 0.4. For even higher shear rates, 
7 > 7**, shear merely reduces the lamellar spacing. The 
significance of these critical values can be discussed in 
relation to the shear rate-dependence of the rheological 
and microscopic properties, which will be discussed else- 
where. 

The parallel sliding layers result from rearrangement of 
the domains of the original zero-rate structures, whereas 
the formation of perpendicular lamellae is accompanied 
by significant breaking-up and remelting of the domains. 
This observation should provide a basis for theoretical 
treatment and provide a subject for further experimental 
study: under shear flow strong enough to overcome the 
segregation tendencies, do block copolymers in a selective 
solvent adopt a perpendicular lamellar structure? 
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